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ABSTRACT 

The Galactic high-mass X-ray binary and jet source (microquasar) 
LSI +61 303 has recently been detected at TeV 7-ray energies by the MAGIC 
telescope. We have applied a time-dependent leptonic jet model to the broadband 
spectral energy distribution and suggested (though not unambiguously detected) 
orbital modulation of the very high energy 7-ray emission of this source. Our 
model takes into account time dependent electron injection and acceleration, 
and the adiabatic and radiative cooling of non-thermal electrons. It includes 
synchrotron, synchrotron self-Compton and external inverse Compton (with seed 
photons from the companion star), as well as 77 absorption of 7-rays by starlight 
photons. The model can successfully reproduce the available multiwavelength 
observational data. Our best fit to the SED indicates that a magnetic field of 
B ~5x 10 3 G at ~ 10 3 R g is required, and electrons need to be accelerated out 
to TeV energies (72 = 10 6 ) with a nonthermal injection spectrum with a spectral 
index of q = 1.7, indicating the operation of acceleration mechanisms beyond the 
standard first-order Fermi mechanism at relativistic or non-relativistic shocks. 
The orbital modulation of the VHE 7-ray emission can be explained solely by 
the geometrical effect of changes in the relative orientation of the stellar com- 
panion with respect to the compact object and jet as it impacts the position and 
depth of the 77 absorption trough. Such a scenario predicts a trend of spectral 
hardening during VHE 7-ray low orbital phases. 
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1. Introduction 

X-ray binaries with relativistic jets, or microquasars, have recently been established as a 
new class of 7-ray emitting sources. Two sources detected by EGRET on board the Compton 
Gamma-Ray Observatory were found to be spatially consistent with the locations of micro- 
quasars, namely LS 5039 and LSI +61 303 (Hartman et al. 1999). Observations of > 250 GeV 
7-rays from LS 5039 with the High Energy Stereoscopic System (HESS) have recently shown 
this high-mass source to be a TeV emitter (Aharonian et al. 2005), also strengthening its ear- 
lier tentative identification with the EGRET source 3EG J1824-1514 (Paredes et al. 2000). 
The second source, LSI +61 303, associated with the COS-B source 2CG 135+01 (Hermsen 
et al. 1977; Gregory & Taylor 1978) and the EGRET source 3EG J0241+6103 (Kniffen et 
al. 1997), was recently observed with the Major Atmospheric Gamma-ray Imaging Chrenkov 
(MAGIC) telescope by Albert et al. (2006). They detected variable 7-ray emission above 100 
GeV over six orbital cycles, suggesting a periodic modulation on the time scale of the orbital 
period of 29.496 d, though further observations are necessary to firmy establish the corre- 
lation between the VHE variability and the orbital period. The strongest detections were 
not at the periastron, suggesting that the modulation might not be related to a modulation 
of the accretion rate (as suggested by, e.g., Romero et al. 2003; Bosch- Ramon & Paredes 
2004b), but rather to geometrical effects. Geometrical effects causing an orbital modulation 
of the high-energy emission include the azimuthal-angle dependence of the 77 absorption 
of high-energy emission by starlight photons (Bottcher & Dermer 2005; Dubus 2006) or of 
the Compton upscattering of starlight photons by relativistic electrons in the microquasar 
jet (Dermer & Bottcher 2006). Dubus (2006) has performed a detailed analysis of the 77 
absorption of VHE photons in the photon field of the stellar companion, taking into account 
the finite size of the star and the eccentricity of the orbit for the case of VHE emission from 
the surface of a putative neutron star associated with LSI +61 303. Note that he used a 
definition of the phase such that phase corresponds to the periastron passage, so that his 
results should be shifted by Aip = 0.23 (see below) when compared to the definition com- 
monly used throughout the literature on this object, including this Letter. He found that the 
77 absorption depth is expected to show a pronounced maximum at phases just before the 
periastron passage, when the compact object is located behind the stellar companion. The 
recent VHE detections, along with the observation of X-ray jet structures in several micro- 
quasars using Chandra and X MM — Newton (Corbel et al. 2002; Tomsick 2002, e.g.,), have 
re-ignited interest in jet models for the high energy emission from microquasars, analogous 
to the commonly adopted models for blazars (for a recent review see, e.g. Bottcher 2002). 

In the leptonic model of microquasars, Very High Energy (VHE) emission might most 
likely originate near the base of the mildly relativistic jet. The soft photons from the com- 
panion star, the accretion disk as well as from jet synchrotron radiation can be Compton 
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upscattered by the ultra-relativistic electrons in the jet. Steady-state leptonic jet models of 
the 7-ray emission from microquasars have been presented, e.g., by Bosch- Ramon & Paredes 
(2004a,b), and (Dermer & Bottcher 2006). A time-dependent, broadband leptonic jet model 
was presented in Gupta et al. (2006), (henceforth, Paper 1) where an analytical solution to 
the electron kinetic equation was presented, restricting the analysis to Compton scattering in 
the Thomson regime. This Letter follows up on the analysis of Paper 1, now incorporating 
a full Klein-Nishina treatment of the Compton scattering as well as the angle dependence of 
the stellar radiation field. We then apply the new model to broadband observations of LSI 
+61 303, iincluding the effect of the orbital modulation on the 77 absorption in the >100GeV 
range. In §2, we present a general outline of the model geometry and the radiative processes 
involved. §3 shows the application to LSI +61 303, and we conclude with a brief summary 
in §4. 



The accretion flow onto the central compact object is ejecting a twin pair of jets, assumed 
to be oriented perpendicular to the orbital plane, which is inclined with respect to the line of 
sight by an angle i. Two intrinsically identical disturbances, containing non-thermal plasma 
(blobs) originate from the central source at the same time, traveling in opposite directions 
along the jet at a constant speed Vj = f3j c. Over a limited range in distance x < x < x ± , 
relativistic electrons are accelerated and injected in the emission region with an exponentially 
cut-off power-law distribution in electron energies (E e = r ym e c 2 ) with low- and high-energy 
cutoffs 71 and 72, respectively, in the co-moving frame: 



The normalization factor Q 1 ^ is approximately related to the power injected into relativistic 
electrons, L in j through 



The injection luminosity L m j and the functional dependence of Eq. (1) are held constant 
between x and 27. The blob's (transverse) radius, R±, scales with distance from the central 
engine as R± = R^i^x/xo) , i.e., a = corresponds to perfect collimation, and a = 1 
describes a conical jet. Following the arguments given in Atoyan & Aharonian (1997), 



2. Model Description 



7i < 7- 



(1) 
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Fig. 1. — Sketch of the geometry of the binary orbit of LSI +61 303, approximating the 
massive (Be) star to be stationary. The labels on the compact-object orbit indicate the 
orbital phase ip as well as the azimuthal angle with respect to the apastron, x- 

we choose a magnetic- field dependence on distance from the central black hole as B(x) = 
B (R ± /R° ± )- 2 = B (x/x )- 2a . 

The geometry of the orbit of the LSI +61 303 system is illustrated in Fig. 1. In order 
to incorporate the orbital modulation of the Compton scattering, we replaced the circular 
orbit of the binary system in Paper 1 by an elliptical one, so that the distance of the star to 
the compact object now varies as 

/ \ ail — e 2 ) 

r(x) = t-t 3 

1 — ecosiv 



where a is the semi-major axis, e is the eccentricity of the orbit, and x is the phase angle 
(which is not linearly related to the orbital phase if), in contrast to the incorrect expression 
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in Romero et al. 2005). The 77 opacity as a function of 7-ray photon energy and orbital 
phase has been calculated following Bottcher & Dermer (2005). 

To treat the time-dependent electron dynamics and the radiation transfer in the emitting 
blob, we solve the continuity equation for the relativistic electrons, 



where, (d^/dt)i oss is the (radiative and adiabatic) energy loss rate for the electrons. The 
electron cooling rates are calculated using subroutines used in the jet radiation transfer code 
of Bottcher, Mause & Schlickeiser (1997) and Bottcher & Bloom (2000). The discretized 
electron continuity equation can be written in the form of a tridiagonal matrix (Chiaberge 
& Ghisellini 1999), and solved using the standard routine of Press et al. (1992). 

Radiation mechanisms included in our simulations are synchrotron emission, Compton 
upscattering of synchrotron photons, namely synchrotron self-Compton (SSC) emission, and 
Compton upscattering of external photons. The dominant external soft photon source is 
the companion star (for details, refer to Paper 1) which is approximated as a point source, 
emitting a blackbody with dimensionless temperature 9* = kT*/m e c 2 . 



The high mass XRB, LSI +61 303, at a distance of 2 kpc (Hutchings & Crampton 
1981; Frail & Hjellming 1991) has been observed from radio (Paredes et al. 1996) to TeV 
7-rays (Albert et al. 2006), providing a wealth of broadband data to test various emission 
models. However, most of the models until now dealt more with the radio, X-ray, or 7- 
ray variability of the source, than with its broadband properties and its very-high-energy 
emission. Recently, however, Chernyakova et al. (2006) reproduced the broadband spectrum 
of the system, including the recent MAGIC detection, using a model of a rotation powered 
pulsar. The (slight) orbital modulation of the X-ray and 7-ray emission in the XMM-Newton 
and INTEGRAL energy ranges was explained based on a variable injection rate of high 
energy electrons. The VHE 7-ray emission was attributed to hadronic processes, initiated 
by high-energy protons in the pulsar wind, and its orbital modulation was not explained in 
that model. In contrast, in this Letter, we are specifically addressing the broadband spectral 
characteristics and the orbital modulation of the VHE 7-ray emission. 





3. 



Application to LSI +61 303 



3.1. 



Parameter Selection 
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Casares et al. (2005) have found a mass function of f(M) 0.011 M Q . Since the 
inclination angle % between the line of sight and the normal to the orbital plane (which, in 
our model is identical to the jet axis) is poorly constrained, and also the mass of the BO V 
type stellar companion (with a dense equatorial wind) is not well known and in the range of 
M* ~ 10 - 15 Mq, a rather wide range of masses for the compact object is consistent with 
the observational constraints, even including a 1.4 M Q neutron star (Hutchings & Crampton 
1981). In our model, we follow the suggestion of Casares et al. (2005) of % — 30° and adopt a 
companion mass of ~ 12M , which would yield M x ~ 2.67 M Q . The orbit has an eccentricity 
of e = 0.72 ± 0.15, the orbital period is ~ 26.496 d (Gregory 2002), and the orbital semi 
major axis is a = 5 x 10 12 cm. MERLIN observations by Massi et al. (2004) imply a jet 
Lorentz factor of Tj et = 1.25. 

In our fitting procedure, we are varying the following parameters: the low and high 
energy cut-off of the electron injection spectrum 71 and 72, respectively, the injection spectral 
index q, the initial distance of the injection zone from the compact object xq, and the initial 
magnetic field Bq. The straight power-law shape of the X-ray spectrum indicates that any 
disk blackbody component from an accretion disk around the compact object should have 
a luminosity of Lp < 10 34 ergs s _1 and be at least ~ 4 orders of magnitude lower than 
the luminosity of the companion star. For such low disk luminosities, external Compton 
scattering of accretion disk photons will be negligible as well. The evolution of the electron 
distribution in the emitting region is followed over a period of ~ 5 days (<C P or b)- Since 
most of the high-energy emission emanates from the system within the first few seconds of 
electron evolution, the orbital phase is taken to be a constant for each run of our model. 
Also, the centre of mass of the system lies close to the massive (12 Mq) stellar companion, 
so that the compact object essentially orbits the stellar companion, which is assumed to be 
stationary. 

For an initial spectral fit, we concentrate on the orbital phase during which MAGIC 
significantly detected the source. Specifically, we chose ip = 0.5, corresponding to a phase 
angle x — —20°. This resulted in a best fit for which all physical parameters are summarized 
in Table 1. 



3.2. Results 

In Fig. 2, we show the computed SEDs of LSI +61 303 using the parameters listed 
in Table 1. The black solid line corresponds to our best fit to the SED at the orbital 
phase of ip = 0.5 (x = —20°). The figure also shows the individual contributions of the 
synchrotron, stellar, and SSC emission, as well as the inverse Compton (EC) scattering of 
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starlight photons. 

The overall agreement of the simulated SED at ip — 0.5 with the observed spectrum is 
very good, given that the multiwavelength data is not simultaneous. Since the radio emission 
most likely originates from regions of the jet further away from the compact object than the 
evolution modeled here, our model underpredicts the VLA data (not included in the figure). 
The NIR data (Strickman et al. 1998) has been de-reddened using N H = 8.4 x 10 21 cm~ 2 
(Taylor et al. 1996). The infrared emission might be dominated by emission from the outer 
parts of the circumstellar disk around the Be star which is not included in our model. 

In order to reproduce the keV - MeV data, a magnetic field of 5 x 10 3 G is required 
at a distance of ~ 1.7 x 10 9 cm from the compact object. In the X-ray band, most of 
the radiation is synchrotron emission from the jet (also see, e.g. Paredes et al. 2006), with 
negligible contributions from SSC or EC. The SSC contribution dominates the spectrum at 
TeV energies. Note that the individual SSC component in Fig. 2 does not include the 77 
absorption feature, which is only applied to the total emanating spectrum. 

After our fit to the ip = 0.5 phase spectrum, we performed a model simulation with 
the same intrinsic jet parameters, but changed the orbital geometry to ip = 0.16 (x = 
70°), corresponding to the MAGIC low state. Since the EC (star) component was rather 
insignificant, this will impact essentially only the 77 absorption. The corresponding model 
fit is illustrated by the solid red line in Fig. 2. In this case, the 77 absorption trough is deeper 
because (a) VHE 7-rays pass by the star more closely than their nearest approach in the ip = 
0.5 geometry, and (b) starlight photons intercept VHE photons at a more favorable angle for 
77 absorption, leading to a lower energy threshold for this process. In this phase, our model 
predicts an integrated flux over 400 GeV of F(E > 400GeV) « 2.2 x 10~ 12 ph cm~ 2 s _1 . This 
is in perfect agreement with the 2a upper limit of F(E > 400GeV) < 3 x 10~ 12 ph cm~ 2 s _1 
given in Albert et al. (2006) and indicated by the maroon arrow in Fig. 2. Consequently, we 
find that the orbital modulation of the VHE emission in LS +61 303 can be explained solely 
by the effect of the azimuthal asymmetry of 77 absorption feature due to starlight photons. 



4. Summary and Conclusions 

We have applied a leptonic jet model to the broadband spectral energy distribution and 
orbital modulation of the VHE 7-ray emission of the microquasar LSI +61 303, taking into 
account time dependent electron injection and acceleration, and the adiabatic and radiative 
cooling of non-thermal electrons. Our model includes synchrotron, synchrotron self-Compton 
and external inverse Compton emission (with seed photons predominantly from the compan- 
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ion star), as well as 77 absorption of 7-rays by starlight photons. Compton scattering is 
treated using the full Klein Nishina cross section and the full angular dependence of the 
stellar radiation field. 

The model can successfully reproduce the available multiwavelength observational data, 
including the most recent MAGIC detection in the TeV range (Albert et al. 2006). Our best 
fit to the SED indicates that a magnetic field of B ~ 5 x 10 3 G at ~ 10 3 R g is required, 
and electrons need to be accelerated out to TeV energies (72 = 10 6 ) with a nonthermal 
injection spectrum with a spectral index of q = 1.7. Such an injection spectrum can not 
be achieved by the first-order Fermi mechanism at relativistic or non-relativistic shocks 
(e.g., Gallant et al. 1999; Achterberg et al. 2001) and therefore suggests that 2 nd order 
Fermi acceleration (Virtanen & Vainio 2005) and/or acceleration at shear boundary layers 
(Ostrowski 2000; Stawarz & Ostrowski 2002; Rieger & Duffy 2004) may play a significant 
role in the acceleration of relativistic particles in microquasar jets. In our model, the X-ray 
and < 1 GeV 7-ray emission is dominated by synchrotron emission from the jet, while SSC 
emission dominates the VHE 7-ray emission. 

The suggested orbital modulation of the VHE 7-ray emission can be explained solely 
by the geometrical effect of changes in the relative orientation of the stellar companion 
with respect to the compact object and jet as it impacts the position and depth of the 77 
absorption trough (Bottcher & Dermer 2005). If this interpretation is correct, a general trend 
of spectral hardening at TeV energies during VHE 7-ray low phases is predicted. Although 
additional effects of a varying accretion rate due to the highly eccentric orbit of the binary 
system are likely to occur, they are not necessary to explain the observed orbital modulation 
of the VHE emission of LSI +61 303. 

We thank the anonymous referee for useful comments and a quick review, and C. 
D. Dermer, G. Dubus, and V. Bosch-Ramon for helpful comments and discussions. This 
work was partially supported by NASA through INTEGRAL GO (Theory) grant award no. 
NNG 05GK59G. 
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Table 1. Relevant parameter choices for our best fit to the broadband spectrum of 

LSI +61 303 



Parameter 


Symbol 


Value 


Distance 


d 


6.17 X 10 21 cm 


Jet inclination angle 


i 


30° 


Bulk Lorcntz factor 




1.25 


Semi Major Axis 


a 


5 X 10 12 cm 


Luminosity of companion star: 




2 X 10 38 ergs s" 1 


Surface temperature of the companion star 


T* 


2.25 X 10 4 K 


Initial blob radius 


Ro 


10 3 Rg 


Jet collimation parameter 


a 


0.3 


Electron injection spectrum, low-energy cutoff 


O'min 


10 


Electron injection spectrum, high-energy cutoff 


7nias 


10 6 


Electron injection spectrum, spectral index 


1 


1.7 


Beginning of electron injection zone 


X 


10 3 Rg 


End of electron injection zone 


XI 


10 5 Rg 


Magnetic field at xq 


Bo 


5 x 10 3 G 


Injection luminosity 




10 35 ergs s _1 


Orbital Period 


P 


26.496 days 


Orbit eccentricity 


e 


0.72 
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v[Hz] 



Fig. 2. — Spectral energy distribution of LSI +61 303. The plotted observational data 
corresponds to if) ~ 0.5 and are taken from: Albert et al. (2006) (MAGIC), Harrison et al. 
(2000) (RXTE PCA + HEXTE), Strickman et al. (1998) (OSSE, NIR), Goldoni & Mereghetti 
(1995) (ROSAT), Leahy (2003) (EGRET), Chernyakova et al. (2006) (XMM, INTEGRAL). 
Also included is the MAGIC 2 o upper limit for phase ip ~ 0.2 (Albert et al. 2006). The solid 
curves show our fit results with 77 absorption calculated for ip ~ 0.5 (black) and ip ~ 0.16 
(red). A trend of spectral hardening at TeV energies during VHE low phases is predicted. 



